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Abstract 

Given a set of positive integers A = {ai, . . . ,a n }, we study the number PA(t) of 
nonnegative integer solutions (mi, . . . , m n ) to J2j=i m j a j = ^ We derive an explicit 
formula for the polynomial part of pa- 



Let A = {cti, . . . , a n } be a set of positive integers with gcd(ai, . . . , a n ) = 1. The 
classical Frobenius problem asks for the largest integer t (the Frobenius number) such 
that 

miai H h m n a n = t 

has no solution in nonnegative integers m\, . . . ,m n . For n = 2, the Frobenius number is 
(ai — l)(a2 — 1) — 1, as is well known, but the problem is extremely difficult for n > 2. 

is to 



(For surveys of the Frobenius problem, see [RjEi].) ^ ne a PP roacxl |BDRt IT| \K\ ISO 

study the restricted partition function the number of nonnegative integer solutions 

(mi, . . . , m n ) to Y^j=i m j a j = where t is a nonnegative integer. The Frobenius number 
is the largest integral zero of PA(t). Note that, in contrast to the Frobenius problem, in 
the definition of pa we do not require a\, . . . , a n to be relatively prime. In the following, 
ai, . . . , a n are arbitrary positive integers. 
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It is clear that Pa{1) is the coefficient of z l in the generating function 

G(z) ' 



(1 - 2 ai ) • • - (1 - Z a ")' 
If we expand G(z) by partial fractions, we see that pA{t) can be written in the form 

A 

where the sum is over all complex numbers A such that \ ai = 1 for some i, and PA,\{t) is 
a polynomial in t. The aim of this paper is to give an explicit formula for PA,i(t), which 
we denote by -Pa(^) and call the polynomial part of Pa{P)- It is easy to see that -Pa(^) is 
a polynomial of degree n — 1. (More generally, the degree of PA,\if) is one less than the 
number of values of i for which X ai = 1.) It is well known jPSl Problem 27] that 

j-n— 1 

(n - 1)! ai ■ ■ -a n 

Our theorem is a refinement of this statement. We note that Israilov derived a more 
complicated formula for Pa(^) in 

Let us define Qyi(t) by Pa(^) = -Pa(^) + From the partial fraction expansion 

above, it is clear that Qa (and hence also pa) is a quasi-polynomial, that is, an expression 
of the form 

c d (t)t d + --- + Cl (t)t + c (t), 

where Co, . . . ,Cd are periodic functions in t. (See, for example, Stanley [EH Section 4.4], 
for more information about quasi-polynomials.) In the special case in which the Oj are 
pairwise relatively prime, each PA,\(t) for A ^ 1 is a constant, and thus Qa(^) is a 
periodic function with average value 0, and this property determines QA(t), and thus 
Pa(^)- Discussions of Qa(£) can be found, for example, in |BDR| H) IK] . 
We define the Bernoulli numbers Bj by 

r.3 



e z - 1 ^ ? 

j>o J 

(so i?o = 1, -£>i = — §, B 2 = |, -B4 = — and S n = if n is odd and greater than 1.} 
Theorem. 



m)! 

m=0 v ' fciH \-k n =m 

i (-i) m 

^•••^^(n-l-m)! 

^ fci!--.Avn! U-llJ '"Im-mlJ ' lJ 

fci+2fc 2 ---+mfc m =m i m \ / \ / 
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where Sj = a\ + • • • + a l n . 

As an intermediate step, we first prove the following more elegant but less explicit 
formula for P A {t)- 

Proposition. P A (t) 2s the constant term in z in 



ze tz 



(1 - e ai2 ) • ■ ■ (1 - e a " 2 ) 
Proof. As noted earlier, is the coefficient of z l in the generating function 

G{z) ' 



(1 - z a ^) •••(!- z a ") 

Hence if we let f(z) = G(z)/z t+1 then p A (t) = Res(/(z), z = 0). As in jBDRj . we use the 
residue theorem to derive a formula for p A (t). Since clearly lim / f(z)dz — 0, 

R ^°°J\z\=R 

PA (t) = - Res(/(z), 3 = 1)-^ Res(/(z), z = A). 

Here the sum is over all nontrivial Oi, . . . , a n th roots of unity A. It is not hard to see 
that Kes(f(z),z = A) may be expressed in the form it A (£)A~* for some polynomial U\(t), 
and thus it follows from our earlier discussion that — Res(/(z), z = \) = P^-i^A - *. In 
particular, 

P A (*) = -Res(/(*),z = l). 
To compute this residue, note that 

Res(f(z),z= 1) = Res(e z f(e z ),z = 0), 

so that 

P A (t) = - Res (- r,* = 0). (4) 

V(l - e aiZ ) • •■ (1 - e a ^) / v ; 

□ 

Proof of the theorem. The coefficient of £ ra_1-m in Pa(£) is by (@J) the coefficient of z~ n+m 
in 

(-l) n " m 1 



(n-l-m)! (1 - e ai2 ) • • ■ (1 - e a 
which is the coefficient of z m in 

f-l) m 



(n — 1 — m)! a\ ■ ■ ■ a n 



B{axz) ■ ■ ■ B(a n z), (5) 



where B(z) = z/(e z — 1), and this implies (J2J). 
To prove (JSJ), we first note that 
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as can easily be proved by differentiating both sides. Then 
B(a lZ ) ■ ■ ■ B{a n z) = exp Vi IV 

-nox P ((-iy'^|). (6) 

Since B 2i+l = for i > 0, (-l)^ 1 ^ = -Bj for j > 1, and © follows from © and ©. 

□ 

Remark. It is possible to avoid the use of complex analysis and give a purely formal power 
series proof of the theorem. We indicate here how this can be done. We work with formal 
Laurent series, which are power series with finitely many negative powers of the variable. 
If F(z) = u i z% is a formal Laurent series (v,i is nonzero for only finitely many negative 
values of i) then the residue of F(z) is resF(2;) = U-i. An elementary fact about formal 
Laurent series is the change of variables formula for residues: If g(z) is a formal power 
series with g(0) = and ^'(O) ^ then 

res.F(,2) = res F(g(z))g'(z). 

(See, for example, Goulden and Jackson |GJt p. 15].) 
By partial fractions, we have 

G(z) = 7 — —, r = + ■■■+ , ° m , + R(Z), 

KJ {I- zl) ■■■{!- z^) 1-z (l-z) m yh 

where R(z) is a rational function of z with denominator not divisible by 1 — z. It follows 
from our earlier discussion that 



E P A(ty = _ii_ + ... + 



Cl 

l. -a 

and thus 



where we take q to be for I > m. Now let U(z) = G(l — z). Then 

1 



U(z) 



(i - (i - z) a ^) • ■ ■ (i - (i - z y™) 

C ± + ... + ^+R(l- Z ), 

z z m 



where R(l— z) has a formal power series expansion (with no negative powers of z), and thus 
q = res z l ^ 1 U(z). Note that this holds for all I > 1, since for / > m, q = res z l ~ 1 U(z) = 0. 
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Then 

U(z) 



1=1 v 7 1=1 v 7 v 



We now apply the change of variables formula with g(z) = 1 — e z and we obtain 

U(l -e z ) 



P A (t) = - res 



e -tz 



— res 



1 - e a i 2 ) • - • (1 - e a ™ z ) 
which is (jlj), and the proof continues as before. 

References 

[BDR] M. Beck, R. Diaz, S. Robins, The Frobenius problem, rational polytopes, and 
Fourier-Dedekind sums, to appear in J. Number Theory. 

[GJ] I. P. Goulden and D. M. Jackson, Combinatorial Enumeration, Wiley, New York, 
1983. 

[I] M.I. Israilov, Numbers of solutions of linear Diophantine equations and their appli- 
cations in the theory of invariant cubature formulas, Sibirsk. Mat. Zh. 22 (1981), no. 2, 
121-136, 237. English translation: Siberian Math. J. 22 (1981), no. 2, 260-273. 

[K] T. Komatsu, The number of solutions of the Diophantine equation of Frobenius, to 
appear in J. Theor. Nombres Bordeaux. 

[PS] G. Polya and G. Szego, Aufgaben und Lehrsdtze aus der Analysis, Springer- Verlag, 
Berlin, 1925. 

[R] J. L. Ramirez Alfonsin, The diophantine Frobenius problem, Report No. 00893, 
Forschungsinstitut fin diskrete Mathematik, Universitat Bonn (2000). 

[Se] E. S. Selmer, On the linear diophantine problem of Frobenius, J. reine angew. Math. 
293/294 (1977), 1-17. 

[SO] S. Sertoz, A. Ozluk, On the number of representations of an integer by a linear form, 
Istanbul Univ. Fen Fat Mat. Berg. 50 (1991), 67-77. 

[St] R. P. Stanley, Enumerative Combinatorics, Vol. 1, Wadsworth & Brooks/Cole, Mon- 
terey, California, 1986. 



THE ELECTRONIC JOURNAL OF COMBINATORICS 8 (2001), #N7 



5 



